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In this paper, the Weiss-Weinstein bound is analyzed in the context of sources localization with a planar array 
of sensors. Both conditional and unconditional source signal models are studied. First, some results are given in the 
multiple sources context without specifying the structure of the steering matrix and of the noise covariance matrix. 
Moreover, the case of an uniform or Gaussian prior are analyzed. Second, these results are applied to the particular 
case of a single source for two kinds of array geometries: a non-uniform linear array (elevation only) and an arbitrary 
planar (azimuth and elevation) array. 
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I. Introduction 

Sources localization problem has been widely investigated in the literature with many applications such as radar, 
sonar, medical imaging, etc. One of the objective is to estimate the direction-of-arrival (DOA) of the sources using 
■ an array of sensors. 

In array processing, lower bounds on the mean square error are usually used as a benchmark to evaluate 
the ultimate performance of an estimator. There exist several lower bounds in the literature. Depending on the 
assumptions about the parameters of interest, there are three main kinds of lower bounds. When the parameters 
are assumed to be deterministic (unknown), the main lower bounds on the (local) mean square error used are the 
well known Cramer-Rao bound [2| and the Barankin bound [3| (more particularly their approximations |4| [5| |6| 
[7 1 [8 1). When the parameters are assumed to be random with a known prior distribution, these lower bounds on 
the global mean square error are called Bayesian bounds Q. Some typical families of Bayesian bounds are the 
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Ziv-Zakai family 02) OH 02 and the Weiss- Weinstein family 03 QH 03 OH- Finally, when the parameter 
vector is made from both deterministic and random parameters, the so-called hybrid bounds have been developed 

H7i ifi8i ifi9i tan. 

Since the DOA estimation is a non-linear problem, the outliers effect can appear and the estimators mean square 
error exhibits three distinct behaviors depending on the number of snapshots and/or on the signal to noise ratio(SNR) 
ll2D . At high SNR and/or for a high number of snapshots, i.e., in the asymptotic region, the outliers effect can 
be neglected and the ultimate performance are described by the (classical/Bayesian/hybrid) Cramer-Rao bound. 
However, when the SNR and/or the number of snapshots decrease, the outliers effect lead to a quick increase of 
the mean square error: this is the so-called threshold effect. In this region, the behavior of the lower bounds are 
not the same. Some bounds, generally called global bounds (Barankin, Ziv-Zakai, Weiss-Weinstein) can predict the 
threshold while the others, called local bounds, like the Cramer-Rao bound or the Bhattacharyya bound cannot. 
Finally, at low SNR and/or at low number of snapshots, i.e., in the no-information region, the deterministic bounds 
exceed the estimator mean square error due to the fact that they do not take into account the parameter support. 
On the contrary, the Bayesian bounds exploit the parameter prior information leading to a "real" lower bound on 
the global mean square error. 

In this paper, we are interested in the Weiss-Weinstein bounds which is known to be one of the tightest Bayesian 
bound with the bounds of the Ziv-Zakai family. We will study the two main source models used in the literature 1 22 1 : 
the unconditional (or stochastic) model where the source signals are assumed to be Gaussian and the conditional (or 
deterministic) model where the source signals are assumed to be deterministic. Surprisingly, in the context of array 
processing, while closed-form expressions of the Ziv-Zakai bound (more precisizely its extension by Bell et. al. 
Il23l ) were proposed around 15 years ago for the unconditional model, the results concerning the Weiss-Weinstein 
bound are, most of the time, only conducted by way of computations. Concerning the unconditional model, in [24|, 
the Weiss-Weinstein bound has been evaluated by way of computations and has been compared to the mean square 
error of the MUSIC algorithm and classical Beamforming using a particular 8x8 element array antenna. In [25], 
the authors have introduced a numerical comparison between the Bayesian Cramer-Rao bound, the Ziv-Zakai bound 
and the Weiss-Weinstein bound for DOA estimation. In ll26l . numerical computations of the Weiss-Weinstein bound 
to optimize sensor positions for non-uniform linear arrays have been presented. Again in the unconditional model 
context, in [27|, by considering the matched-field estimation problem, the authors have derived a semi closed-form 
expression of a simplified version of the Weiss-Weinstein bound for the DOA estimation. Indeed, the integration 
over the prior probability density function was not performed. The conditional model (with known waveforms) is 
studied only in |28], where a closed-form expression of the WWB is given in the simple case of spectral analysis 
and in [1J which is a simplified version of the bound. 

While the primary goal of this paper is to give closed-form expressions of the Weiss-Weinstein bound for the DOA 
estimation of a single source with an arbitrary planar array of sensors, under both conditional and unconditional 
source signal models, we also provide partial closed-form expressions of the bound which could be useful for 
other problems. First, we study the general Gaussian observation model with parameterized mean or parameterized 
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covariance matrix. Indeed, one of the success of the Cramer-Rao is that, for this observation model, a closed-form 
expression of the Fisher information matrix is available: this is the so-called Slepian-Bang formula [29 1 . Such kind 
of formulas have been less investigated in the context of bounds tighter than the Cramer-rao bound. Second, some 
results are given in the multiple sources context without specifying the structure of the steering matrix and of the 
noise covariance matrix. Finally, these results are applied to the particular case of a single source for two kinds 
of array geometries: the non-uniform linear array (elevation only) and the planar (azimuth and elevation) array. 
Consequently, the aim of this paper is also to provide a textbook of formulas which could be applied in other fields. 
The Weiss-Weinstein bound is known to depend on parameters called test points and other parameters generally 
denoted Si. One particularity of this paper in comparison with the previous works on the Weiss-Weinstein bound 
is that we do not use the assumption Si — 1/2, Vi. 

This paper is organized as follows. Section [EI] is devoted to the array processing observation model which will be 
used in the paper. In Section [TTTJ a short background on the Weiss-Weinstein bound is presented and two general 
closed-form expressions which will be the cornerstone for our array processing problems are derived. In Section 
ITVl we apply these general results to the array processing problem without specifying the structure of the steering 
matrix. In Section [V] we study the particular case of the non-uniform linear array and of the planar array for which 
we provide both closed-form expressions of the bound. Some simulation results are proposed in Section [VT] Finally, 
Section |yil] gives our conclusions. 

II. Problem setup 

In this section, the general observation model generally used in array signal processing is presented as well as the 
first different assumptions used in the remain of the paper. Particularly, the so-called conditional and unconditional 
source models are emphasized. 

A. Observations model 

We consider the classical scenario of an array with M sensors which receives N complex bandpass signals 
s (t) — [si (t) S2 (t) ■ ■ ■ sn (t)] T ■ The output of the array is a M x 1 complex vector y (t) which can be modelled 
as follows (see, e.g., Il30l or [22 1) 

y(t) = A(e)s(t) + n(t), t = l,...,T, (1) 

where T is the number of snapshots, where 9 = [6± 62 ■ ■ ■ q ] T is an unknown parameter vector of interest, where 
A (8) is the so-called M x N steering matrix of the array response to the sources, and where the M X 1 random 
vector n (t) is an additive noise. 

'Note that one source can be described by several parameters. Consequently, q > N in general. 
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B. Assumptions 



• The unknown parameters of interest are assumed to be random with an a priori probability density function 
p(0i) , i = 1, ... ,q. These random parameters are assumed to be statistically independent such that the a 



priori joint probability density function is p (9) = JJp (^*)- We also assume that the parameter space, denoted 

i=l 

Q, is a connected subset of R 9 (see (3T)). 
• The noise vector is assumed to be complex Gaussian, statistically independent of the parameters, i.i.d., circular, 



restrictive in Section [V] where it will be assumed that R n = cr^I. In any case, R n is assumed to be a full 
rank matrix. 

• The steering matrix A (9) is assumed such that the observation model is identifiable. From Section [III] to 
Section [TV] the structure of A (9) is not specified in order to obtain the more general results. 

• Concerning the source signals, two kinds of models have been investigated in the literature (see, e.g., (321 or 
ll22l ') and will be alternatively used in this paper. 

- A4\: Unconditional or stochastic model: s(t) is assumed to be a complex circular random vector, i.i.d., sta- 
tistically independent of the noise, Gaussian with zero-mean and known covariance matrix E [s (t) s H (f )] = 
R s . Note that concerning the previous results on the Cramer-Rao bound available in the literature |32|, 
the covariance matrix R s is assumed to be unknown. In this paper, we have made the simpler assumption 
that the covariance matrix R s is known. These assumptions have already been used for the calculation of 
bounds more complex than the Cramer-Rao bound (see, e.g., [27|, [33 1, |34|). 

- M.2'- Conditional or deterministic model: Vt, s(t) is assumed to be deterministic known. Note that, under 
the conditional model assumption, the signal waveforms can be assumed either unknown or known. While 
the conditional observation model with unknown waveforms seems more challenging, the conditional 
model with known waveforms signals which will be used in this paper can be found in several applications 
such as in mobile telecommunication and radar (see e.g. 11351 . (361 . (37], (38], and [39 1). 

C. Likelihood of the observations 

Let R y = E [y (t) y H (f )] be the covariance matrix of the observation vector y (t) . According to the aforemen- 
tioned assumptions, it is easy to see that under M.\, the observations y (t) are distributed as a complex circular 
Gaussian random vector with zero mean and covariance matrix R y (9) = A(9)H S A H (9) + R„ while under M.2, 
the observations y (t) are distributed as a complex circular Gaussian random vector with mean A(9)s (t) and 
covariance matrix R y = R n . Moreover, in both case the observations are i.i.d.. 

Therefore, the likelihood, p (Y; 9) , of the full observations matrix Y = [y (1) y (2) ... y (T)] under Ai\ is 
given by 



with zero mean and known covariance matrix E [n (t) n 



H 



(t)l = R n . This assumption will be made more 




(2) 
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where R y (0) = A(9)~R S A H (0) + R n and the likelihood under M2 is given by 

p (Y; 0) = 1 j, cxp (y (*) - A (0) s R n 1 (y (t) - A (0) s (*)) j . (3) 

|Rn| V t=l / 

III. Weiss-Weinstein bound: Generalities 

In this Section, we first remind to the reader the structure of the Weiss-Weinstein bound on the mean square 
error and the assumptions used to compute this bound. Second, a general result about the Gaussian observation 
model with parameterized mean or parameterized covariance matrix, which, to the best of our knowledge, does not 
appear in the literature is presented. This result will be useful to study both the unconditional model M.\ and the 
conditional model M.2 in the next Section. 



A. Background 

The Weiss-Weinstein bound for a q x 1 real parameter vector is a q x q matrix denoted WWB and is given 
as follows Il40l 

WWB = HG _1 H T , (4) 

where the q x q matrix H = [hi I12 . . . h 9 ] contains the so-called test-points hi, i = 1, . . . , q such that + hi 6 6 
Vhj. The k, I — element of the q x q matrix G is given by 

E [(L°* (Y; + h fc , 0) - L 1 ~ Sk (Y; - h fc , 6)) (Y; 6 + h t , 8) - L 1 ^ 8 ' (Y; 6 - h t , 0))] 
1 h - 1 E[L«*CY;0 + h fc ,0)]E[L«»(Y;0 + h,,0)] ' () 

where the expectations are taken over the joint probability density function p (Y, 6) and where the function 
L (Y; + h.,,6) is defined by L (Y; 8 + h u 8) = ' ■ The elements s t are such that s t G [0, 1], % = 1, . . . , g. 

Note that we have the following order relation 1401 

(») ' ' 



0-0) [8-0 



y WWB, (6) 



Cov [0j = E 

where A >z B means that the matrix A — B is a semi-positive definite matrix and where Cov [0j is the global 
(the expectation is taken over the joint pdf p (Y, 0)) mean square error of any estimator of the parameter vector 
0. Finally, in order to obtain a tight bound, one has to maximize WWB over the test-points h^ and Si i = 1, . . . , q. 
Note that this maximization can be done by using the trace of HG _1 H T or with respect to the Loewner partial 
ordering 0411 . In this paper we will use the trace of HG _1 H T which is enough to obtain tight results. 

B. A general result on the Weiss-Weinstein bound and its application to the Gaussian observation models 

An analytical result on the Weiss-Weinstein bound which will be useful in the following derivations and which 
could be useful for other problems is derived in this part. Note that this result is independent of the parameter 
vector size q and of the considered observation model. 
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Let us denote the observation space. By rewriting the elements of matrix G (see Eqn. (f5]l) involved in the 
Weiss-Weinstein bound, one obtains for the numerator denoted -W{G} fc , > 

N {G}ki = E [(L Sk (Y; 9 + h k) 0) - L^ Sk (Y; 9 - h k ,0)) (i Sl (Y; 9 + h,, 9) - L 1 ^' (Y; 9 - hi,0))] 

p s * (y, e + h fc )p« (y, + ^\ dYde+ f f p'- Sk (Y,9-h k )p^ Q^g-h,) 



e Jo 



e Jn 



p-*-«(Y,e) y e y p««-»(Y,e) 



and for the denominator denoted D{ G } fc ; , 



£ {G}fc ; = E [L s <* (Y; + h fe , 9)] E [L*< (Y; 9 + hi, 6)] 

pSk ^ 9 + ^ d Y d 9[ K ( ^* + h 'W (8) 



, eJn ^-i(Y,0) -'""J e J Q p*i-i(Y,6) 
Let us now define a function 77 (a, /3, u, v) as 

1 r \ f f P° (Y,9 + u)pP(Y,9 + v) 
77 (a, P , u, v) = J & ^ pa+ ^ [Yj6) dYdB, (9) 

where (a, f3) e [0, l] 2 and where (u, v) are two qxl vectors such that 9 + u 6 6 and 9 + v 6 9. By identification, 
it is easy to see that 

m kjl = 

77 (Sfc, Si, hfc, h;) + 7] (1 - s fc , 1 - si, -h fc , -h;) - T) (Sk, 1 - Sj, hfc, -h;) - 77 (1 - Sfc, -hfc, h;) 

7? (s fc ,o,hfe, 0)77(0,^,0, hi) 

Note that we choose the arbitrary notation -D{G} fc , = ?7 ( s fci 0, hfc, 0) 77 (0, si, 0, h;) for the denominator. The 
notation £>{ G } fc , = v( s k, 1, hfe, 0) 77 (1, s t , 0, h;) or, even, £>{ G } fc , = («fe, 0, h fc , v) 77 (0, s t , u, h ; ) will lead to the 
same result. 

With Eqn. dlOt . it is clear that the knowledge of 77 (a, /3, u, v) for a particular problem leads to the Weiss-Weinstein 
bound (without the maximization procedure over the test-points and over the parameters Sj). Surprisingly, this simple 
expression is given in l40l only for s$ = |, Vi and not for the general case. 

Let us now detail this function 77 (a, f3, u, v). The function 77 (a, (3, u, v) can be rewritten as 

f p a (0 + u)pf (g + v) /■ p° (Y;0 + u)^(Y;0 + v) ^ Q 

77 ft V) = L (0) L P^-HY;9) dYM 

, , „ n P q (0 + u)t/(0 + v) 
Q % («• A U > v ) ^ ^J-i (g) (ID 



where we define 



• fll / (Y;g + u)/ (Y;0 + v) ^ 

77 e (a,/3,u,v,0)= / p«+/9-i(Y;g) dY. (12) 



Our aim is to give the most general result. Consequently, we will focus only on fj e (a, f3, u, v) since the a priori 
probability density function depends on the considered problem. 

An important remark pointed out in OTI is that the integration for the parameter space is with respect to 
the region {9 : p (9) > 0} . However, since the functions being integrated are p (9) , p (9 + u) , and p (9 + v) , 
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then the actual region of integration (where all the functions are positive) is the intersection of three regions, 
{9 : p (9) > 0} n {9 : p (9 + u) > 0} H {9 : p (9 + v) > 0} . Note that, in order to simplify the notation we only 
use throughout this paper but this remark will be useful and explictely specified in Section IIV-BI 

1) Gaussian observation model with parameterized covariance matrix: One calls (circular, i.i.d.) Gaussian 
observation model with parameterized covariance matrix, a model such that the observations y (t) ~ CM (0, R y (9)) 
where 9 are the parameters of interest. Note that Aii is a special case of this model since the parameters of 
interest appear only in the covariance matrix of the observations which has the following particular structure 
R y (9) = A(9)H S A H (9) + R n . The closed-form expression of i) g (a, 0, u, v) is given by: 

f, e (a,0,u,v) = - — ggUl - T . (13) 

|Ry(0 + u)f Q |Ry(0 + v)| T/3 \aH y \9 + u) + /3IL y 1 {9 + v)-(a + 0- l)Ry \9)\ 

The proof is given in Appendix lAl Note that, similar expressions are given in |23l (Eqn. (B.15)) and 11421 (p. 67, 
Eqn. (52)) for the particular case where a = s and = 1 — s. 

2) Gaussian observation model with parameterized mean: One calls (circular, i.i.d.) Gaussian observation model 
with parameterized mean, a model such that the observations y (t) ~ CN (f (9) ,R y ) where 9 are the parameters 
of interest. Note that M.2 is a special case of this model since the parameters of interest appear only in the mean of 
the observations which has the following particular structure i t (9) = A(9)s (t) (and R y = R n ). The closed-form 
expression of rj (a, (3, u, v) is given in this case by 

T 

In r, e P, u, v) = - J^a (1 - a) ff (9 + u) R y H t {9 + u) +/? (1 - 0) ff {9 + v) R y H t (9 + v) 

t=i 

+ (l-a-p)(a + p) t t H (9) Ry-H (9) - 2 Re {a/3ff (9 + u) R y J f f (9 + v) 

+a (1 - a - P) f t H (9 + u) Ry^ft (O) + 0(l-a-0) ff (6» + v) R y (0)} , (14) 

or equivalently by 

T 

kn7 e (a, /3, u, v) - - J^a (1 - a - 0) ||R y ^ (f t (0 + u) - f t (0))|| 2 + a/3 ||R y l ' 2 (f t (0 + u) - f t (0 + v)) 

t=i 

+/3 (1 - a - j8) ||R y 1/2 (f t (0 + v) - ^ (0))|| 2 . (15) 
The details are given in Appendix [B] 

IV. General application to array processing 

In the previous Section, it has been shown that the Weiss-Weinstein bound computation (or, at least, the matrix G 
computation) is reduced to the knowledge of the function r\ (a, 0, u, v) given by Eqn. (0. As one can see in Eqn. 
( flOl l, the elements of the matrix G depend on r\ (a, 0, u, v) for particular values of a, 0, u, and v. Consequently, 
the goal of this Section is to detail these particular functions for our model given by Eqn. (Q]). Since Eqn. (0 can 
be decomposed into a deterministic part (in the sense where r\ e (a,0,u,v) (see Eqn. ( fl2b ) only depends on the 
likelihood function) and a Bayesian part (when we have to integrate rj g (a, 0, u, v) over the a priori probability 
density function of the parameters), we will first focus on the particular functions rj (a, 0, u, v) by using the results 
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of the previous Section on the Gaussian observation model with parameterized mean or covariance matrix. Second, 
we will detail the passage from f\ e (a, (3, u, v) to 77 (a, f3, u, v) in the particular case where p{9i) is a uniform 
probability density function Vi. Another result will also be given in the case of a Gaussian prior. 



A. Analysis of r] e [a, /3, u, v) 

We will now detail the particular functions r\ e (a, f3, u, v) involved in the different elements of {G} k t , k, I 6 
{1, q} 2 for both models Mi and M.2- 

1) Unconditional observation model M.\: Under the unconditional model Mi, by using Eqn. (1131 1. one obtains 
straightforwardly the functions ij g (a,(3,u, v) involved in the elements {G} k ; = {G} ; k 

'is s h h ) — |R y (6>)| T 0fc+ a '- 1 ) 

Mi s k) 1 - ai) -h k , -ho = iR y wr(;-^-)|R y ( e -h fc) r(^- i )|R y ( 9 -h i) r(- : - 1 ) 



- s fc ,s;,-h fe ,h ; ) 
57e(s fc , 0,h fe ,0) 
%(0,Si,0,hi) : 



|R y (6+h fc )| Ts fc |s fc R y - 1 (e+h fc )+(l-s i )Ry ^O-hO-Csfc-sORy 1 W| ' 

|R y (e)r( 3 i- 3 fc)|R y (e-h fc )r( 3 fc- 1 ) 

|R y (e+h ^ )^ 3 '|(l-s fc )R y 1 (e-h fc )+s i R y 1 (e+h ^ )-(s ^ -s fc )R y 1 (e)|' J, ' 



(16) 



|R y (e+h fc )| Ts fc| Sfc R y - 1 (e+h fc )-(s fc -i)R y 1 (e)| T ' 
IRrWr^- 1 ) 



|R y (e+h ! )r s '|s ! R y 1 (e+h i )-( Si -i)R y 1 (e)| 
The diagonal elements of G are obtained by letting k — ( in the above equations. 

2) Conditional observation model M.2'- Under the conditional model M.%, by using Eqn. (O with f t (0) = 
A (0) s (t) and R y = R n one obtains straightforwardly the functions rj (a, (3, u, v) involved in the elements 

\nij g (s k ,si,h k ,h.i) = s k (s k + sj - 1) Ce (hfc,0) + sj (s fc + s; - 1) Ce (h;,0) - s fc SiCe (hfe,h ; ) , 
Info (1 - s fc , 1 - s h -hfc, -h t ) = (s k - 1) (s k + st - 1) C e (-hfc, 0) + (a, - 1) (s fc + s ; - 1) Ce (-hi, 0) 

-(l- S fc)(l- S; )C e (-h fe ,-h ; ), 
Info (sfc,l - s;, hfc, -hi) = s k (sfc - sj) Ce (h fe ,0) + (1 - s t ) (s k - s ; ) Ce (-h;,0) + s k (sj - 1) Ce (h fc , -h,) , 
Info (1 - Sfc, s;, -hfc, h ; ) = (s fc - 1) (s k - s;) Ce (-hfc,0) + s; (s ; - s fc ) Ce (h;,0) + (s k - 1) si( e (-h fc ,h;) , 
Info (sfc,0, hfc,0) = s k (s k - 1) Ce (h*,0) , 
Info (0,s/,0,h/) = s; (s; - 1) Ce (h/,0) , 



where we define 



Ce (m, P) = || R n 1/2 (A(6 + n)-A(0 + p)) s (t) 



t=i 



(17) 



(18) 



The diagonal elements of G are obtained by letting k — I'm the above equations. Note that, since we are working 
on matrix G, all the previously proposed results are made whatever the number of test-points. 
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B. Analysis of rj (a, /3, u, v) with a uniform prior 

Of course, the analysis of ij (a, (3, u, v) given by Eqn. ( fTTT i can only be conducted by specifying the a priori 
probability density functions of the parameters. Consequently, the results provided here are very specific. However, 
note that, in general, this aspect is less emphasized in the literature where most of the authors give results without 
specifying the prior probability density functions and compute the rest of the bound numerically (see e.g., ||27 | [25 1 

ma). 

We assume that all the parameters di have a uniform prior distribution over the interval [di, bj\ and are statistically 
independent. We will also assume one test-point per parameter, otherwise there is no possibility to obtain (pseudo) 
closed-form expressions. Consequently, the matrix H is such that 

H = Diag{[h 1 h 2 ---h q }), (19) 

and the vector hi, i = 1, . . . , q, takes the value hi at the i th row and zero elsewhere. So, in this analysis, the vector 
u takes the value Uj at the i th row and zero elsewhere and the vector v takes the value Vj at the j th row and zero 
elsewhere (of course, we can have i = j). Under these assumptions, r/ (a, /3, u, v) can be rewritteij^ for i ^ j 

i r \ f < ( r \ p a (g«±"<)pg (gj+M^ (gj) rr (a \ ah 

^(a^,u,v) = ^ 9 (a,ftu,v) ^-1^)^-1^.) 11 P^ e ^ d9 

! II %(«,ftu,v),iMy(e/{9„y), (20) 



'i 



fc=l 



[a,, bi — Ui] if Ui > 0, [a,j, bj — Vj] if Vj > 0, 

where Qi = ^ and Qj ■ — < . For i = j, one can have v = ±u, 

[dj - Mi, &»] if ^ < 0, [oj- - Uj, bj] if < 0, 

then one obtains 

i] {a, 0, u, v = ±u) = / r/ e (a, /3, u, v) a+/ = ) „ 1 t [ [P (0k) dd 

Je P \ 1) k=1 



f f f le (a,p,u,v = ±u)d6 i d(e/{9 i }). (21) 

j _[ (Ofe - a k ) 



fe=i 



[di + Ui,bi — if Ui > 0, 

In the last equation, if v = u, then Si = ^ , while, if v = u, then Qi = 

[di - Ui, bi + Ui] if Ui < 0, 

[di, bi - Ui] if Ui > 0, 
[di - Ui,bi] if Ui < 0, 

Depending on the structure of rj (a, (3, u, v) , r\ (a, (3, u, v) has to be computed numerically or a closed-form 
expression can be found. 

2 In this case, one has to have a particular attention to the integration domain as mentionned in Section IIII-BI It will not be the case for the 
Gaussian prior since the support is M. 



November 29, 2012 



DRAFT 



10 



Another particular case which appears sometimes is when the function rj e (a, ft u, v) does not depend on 9 
(see, ESI Il9l IT21 fl23l fl25l flZgl OTJ fl33l and Section [V] of this paper). In this case, rj g (a, ft u, v) is denoted 
77 (a, ft u, v) and one obtains from Eqn. d20i > 



77 (a, ft u, v) = — 

Y[ ( b k - a k) 

k=i 



( • 

n 

k=l 







/ d9 k 


/ dBi f 






) 




T-^Ka,/?,u,v), 



and from Eqn. ( f2TT > 



and 



(h - Oi) (6j - Oj) 



/ « s {bi - <H - \ui\) , , 
V (<*, ft u, v = u) = 17 (a, ft u, v) , 

(bi - di) 



77 (a, ft u, v = -u) = r 77 (a, ft u, v) . 

(0i - Oi j 



(22) 



(23) 



(24) 



C. Analysis 0/77 (a, ft u, v) w'f/i a Gaussian prior 

Finally, one can mention that if the prior is now assumed to be Gaussian, i.e., 0, ~ Af erf) Vz and i) g (a, ft u, v) 
does not depend on 9 one obtains after a straightforward calculation 

P a (6i + Ui) 



r?(a,ftu,v) = J7(a,ftu,v) 







= 77(a,ftu,v)exp -- - 



(25) 



and 



77 (a, ft u, v = u) = 77 (a, ft u, v) / 



p"+0 (gH-uQ 
p«+/s-i (0,) 



77 (a, ft u, v) exp - 



(a + ft (1 - a - /3) 
2rf 



77 (a, ft u, v = -u) = 77 (a, ft u, v) 



P"+' 3 - 1 (0i) 
(a + g - a 2 - p 2 + 2a/?) M 2 
2^? 



77 (a, ft u, v) exp 



(26) 



(27) 



V. Specific applications to array processing: DOA estimation 

We now consider the application of the Weiss-Weinstein bound in the particular context of source localization. 
Indeed, until now, the structure of the steering matrix A (9) for a particular problem has not been used in the 
proposed (semi) closed-form expressions. Consequently, these previous results can be applied to a large class of 
estimation problems such as far-field and near-field sources localization, passive localization with polarized array 
of sensors, or radar (known waveforms). 
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Here, we want to focus on the direction-of-arrival estimation of a single source in the far-field area with narrow- 
band signal. In this case, the steering matrix A (9) becomes a steering vector denoted a(0) (except for one 
preliminary result concerning the conditional model which will be given whatever the number of sources in Section 
|V-A21 >. The structure of this vector will be specified by the analysis of two kinds of array geometry: the non-uniform 
linear array from which only one angle-of-arrival can be estimated (6 becomes a scalar) and the arbitrary planar 
array from which both azimuth and elevation can be estimated (0 becomes a 2 x 1 vector). In any cases, the array 
always consists of M identical, omnidirectional sensors. Both model Aii and A4 2 will t> e considered and the noise 
will be assumed spatially uncorrelated: R n = rr^I. Since we focus on the single source scenario, the variance of 
the source signal s (t) is denoted a 2 s for the model M.\. 

The general structure of the i th element of the steering vector is as follows 



where represents the parameter vector, where A denotes the wavelength, and where r$ denotes the coordinate of 
the i th sensor position with respect to a given referential. In the following, will be a scalar or a 2 x 1 vector 
depending on the context (linear array or planar array). 

A. Preliminary results 

Since our analysis is now reduced to the single source case, we give here some other closed-form expressions 
which will be useful when we will detail the specific linear and planar arrays. 

1) Unconditional observation model M.\: In order to detail the set of functions r] e given by Eqn. ( fTo*b . one 
has to find closed-form expressions of the determinant |R y (0 + u)| and of determinants having the following 
structure: |miR~ 1 {0\) + m2R y 1 (^2 ) | with mi + m>2 = 1 or |miRy 1 (9i) + W2Ry 1 (02) + "^Ry 1 ( ^3 ) | with 
mi + rri2 + 7713 = 1. Under M.\, the observation covariance matrix is now given by 




(28) 



Ry(0) 



(T?a(0)a H (0) + a 2 n I M - 



(29) 



Concerning the calculation of |R y (0 + u)|, it is easy to find 




(30) 



Moreover, after calculation detailed in Appendix ICl one obtains for the other determinants 




1 




(31) 
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and 



m 1 R y 1 (0 1 ) + m 2 R y 1 (6/ 2 )+m 3 R y 1 (6» 3 )| = 
if 3 i 3 3 

pncr ( 1 - 5>*P* IK^f - ^EE™^™^' (||a ff (0 fc )a(^)H 2 - l|a(e fc )|| 2 ||a(0* 



/ 3 \ / 3 333 

(n m *V* nil a (^)H 2 -5EE E ||a ff (0 fc )a(0 fe O|ri|a(0 fe ")l| 2 

\ k' ^kk" ^k' ^k 

-a H (0 3 )a(02)a H (0i)a(0 3 )a H (0 2 )a(0 1 ) + a" (0 3 )a(0 1 )a ff (0 1 )a(0 2 )a ff (0 2 )a(0 3 ))) , 



where 



-, fc=l,2,3. 



(32) 



(33) 



<^||a(e fc )ir + ^ 

2) Conditional observation model 7W 2 : Note that the results proposed here are in the context of any number of 
sources. Under the conditional model, the set of functions r) e given by Eqn. ( fl7l i is linked to the function C, e (p, p) 
given by Eqn. d 1 8t > . In this analysis, the vector p takes the value p, i at the i fh row and zero elsewhere and the 
vector p takes the value pj at the j th row and zero elsewhere (of course, one can has i = j). In Appendix iDl the 
calculation of the following closed-form expressions for ( g (p, p) are detailed. 

• If (m — l)p + 1 < i, j < mp, where p denotes the number of parameters per source, then, we have 



T At At 

Ce (m, p) = E "I 8 (*)U a EE { R n% ex p ( j't ( r J - r D 

t=i »=ij=i 



.2tt 

T 



.2tt 

T 



27T 



Otherwise, if (m — l)p + 1 < i < mp and (n — l)p + 1 < j < rip, then we have 
Ce (Mj p) = 



At At 



i=lj=l 
At At 



E Hi 8 WmW 2 EE l R n 1 } tJ eX P ( 3~C ( r J ~ r D ™ ) CX P ( t* m ) ex P ( Jj^M 



A 



2tt 



A 



.2tt 



E IK* (*)U a EE {R, 1 },,, ex p ( r J - r H »») -p (-^Pn) ex p 

{=1 i=lj=l ^ / V / 

2Reff^{s(t)};{ S (t)}„ 
\t=i 

M M ( 2tt A 

x EE { R »% ex P ( J-r ( r J e » " r ^™) ) ex P 

i=l j = l ^ ' 



A 

.2?r 
T 



. . A 
i=ij= 

In particular, if one assumes R„ = cr^I, then, several simplifications can be done: 
• If (m — l)p + 1 < i, j < mp, then 



J-r 1 "* Mm ex P J-r r 7 P 



(35) 



1 M 

(O (p,p) = -2-E 



27T y 



27T rp 



2 T 



Eii^wi 



2 

mil ' 



(36) 



t=i 
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Fig. 1. 3D source localization using a planar array antenna. 

where we note that the function £ e (/x, p) does not depend on the parameter 6. 
• Otherwise, if (to — l)p + 1 < i < mp and (n — 1) p + 1 < j < rtp, then 



Ce O, P) 



1 M 

^E 



exp 



27T 

-J— r- p r , 



2 T M 

E"{ s WU 2 + -rE 



exp 



2tt rp 
~j — r i Pr 



2 T 



-2 Re 



^E exp (j-J- r f ( 6 » ~ e m) ) j exp ^-jyrf/2 m ^ exp ^jyrfp n ^ E{ s WO s Wj»^ (37) 

It is clear that the proposed above formulas for both the unconditional and the conditional models can be applied 
to any kind of array geometry and whatever the number of sources. However, they generally depend on the parameter 
vector 6. This means that, in general, the calculation of the set of functions 77 will have to be performed numerically 
(except if one is able to find a closed-form expression of Eqn. (fTTTl). In the following we present a kind of array 
geometry where, fortunately, the set of functions rj g will not depend on 6 leading to a straightforward calculation 
of the bound. 



B. 3D Source localization with a planar array 

We first consider the problem of DOA estimation of a single narrow band source in the far field area by using 
an arbitrary planar array. In fact, we start by this general setting because the non-uniform linear array is clearly a 
particular case of this array. Without loss of generality, we assume that the sensors of this array lay on the xOy plan 
with Cartesian coordinates (see Fig. [T). Therefore, the vector contains the coordinate of the i th sensor position 
with respect to this referential, i.e., = [d Xi d yi ] T , i = 1, . . . , M. From 



a(0) 



ex P (j^(d Xl u + d Vi v)j ... exp (j^-(d XM u + d yM v)^j 



the steering vector is given by 

n T 



(38) 



where, as in (23], the parameter vector of interest is 9 = [u v] T where 



(39) 



u = sin ip cos (f>, 
v = sin ip sin (p, 

and where <p and </> represent the elevation and azimuth angles of the source, respectively. The parameters space 
is such that u G [—1,1] and v e [—1,1]. Therefore, we assume that they both follow a uniform distribution over 
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[—1,1]. Note that from a physical point of view, it should be more tempting to choose a uniform prior for ip 
and (f>. This will lead to a probability density functions for u and v not uniform. To the best of our knowledge, 
this assumption has only been used in the context of lower bounds in [25 1. Unfortunately, such prior leads to an 
untractable expression of the bound (see Eqn. (21) of [25]). Consequently, other authors have generally not specified 
the prior leading to semi closed-form expressions of bounds (i.e. that it remains a numerical integration to perform 
over the parameters) ll25l ||431 11271 . On the other hand, in order to obtain a closed-form expression, authors have 
generally used a simplified assumption, i.e. a uniform prior directly on u and v (see, for example, [26] [44|). In 
this paper, we have followed the same way by expecting a slight modification of performance with respect to a 
more physical model and in order to be able to get closed-form expressions of the bound. 
We choose the matrix of test points such that 

r h u 







(40) 



Then, we have: 9+h u = [u + h u v] T and 9+h v = [u v + h v ] T . Moreover, we now have two elements 
Si G [0, 1] , i = 1, 2 for which we will prefer the notation s u and s v , respectively. 

1) Unconditional observation model M.\: Under M.\, let us set Usnb, — -2 tttA The closed-form 

{G} u 
{Gfvu {Gjv 



expressions of the elements of matrix G = 
/ 



are given by (see Appendix [E] for the proof): 



(l - % i ) 1 + 2*«(1 - 2s u )Usnr M 



+ ^) 1 + 2(1 - s u )(2s v - 1)Usnr M 



M 

fc=i 
2 



\ 



M 

E ex P (-ifiA) 
fc=i 



{Gh 



2(1- \hu\) l + s„(l -s u )Usnr M 



M 

E ex P {-3^-d Xk h u ) 



1_^M 1 + Su (i _ Su )U t 



SNR 



M 2 



M 



E en>{-j-£d Xk h u 



k=l 



-2T 



(41) 



(i + 2s v (1-2s v )Usnr (m 



M 

E 

fc=i 



+ 1 - \Hf 1 + 2(1 - b„)(2*„ - 1)Usnr M 



E e*P{-3^dy k h v ) 



{G} vv 



2(l-|ftt,|) l+s v (l- s v )Usnr M 



E exP {-j^r d y k h A 



( l - l -¥) 2 U + s v (i-s v )u SN R Up 



E 

fc=l 



E exp (-Jx^Al 



-2T 



(42) 
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E exp(-j2f {d Xk h u -d yk h v ))\ 

k = l 

/||M 

+s u (l-s u -s„) ( || £ exp(-j^d Xk h u )\ 

- s u - 



e ° x p (-i^f^i/fc' 1 ") 
i j==i. 



M 

- M 

- M- 



tt'2 



-S U S V (1 - 8„ - S„) SW | 



V 



£ exp (.-E^W f; cxp £ exp /. ^».-^>..) 

fc = l v ' fc=l v ' k = l \ 



M . 2-Kd y h v \M ( .Z-*d x h u \ ™ ( 

+ E CX P ( -J — T — ) E ox P ( 3 — r — ) 22 ex P - 

k=l v ' k=l v ' fc=l V 



.2w(d X . h U -dy. hy) 



-M 
— M 

1 - Usis 



f2 exp (-j^f dyfe'ii.) - M E cxp (~i^f ^ fe 'i tl ) 



E exp(-.j2f {d x .h u - d y h v ))\\ + M 3 
k=i 



(1 - s„)(l - s„) E exp(j2f (d^/i,, - d yfc h„)) - M- 
\ II fc=i 

§1 M II 2 
E exp (jSf - M< 

I M II 2 
E exp(j2f d Sfc /i„) - Af 2 
I fc=l 



(1 - s u )(l - s v )(s u + 



1) 



S.Yj.'" 



E cx p ( J — t — ) E ex P ( -3 — r — ) E ex P 3 is — i — - — I 

k=i v ' k=i v ' k=i V / 

+ E exp {-j ^ ) Yl exp (j ^ J E exp I -j ^-^ ^ I 

fc=l fc=l fc=l \ / 



-Af 
-Af 



M || 2 II M | 

E exp (-j^dyfcfcv) - M E CX P (-3^-d Xk h u )\ 

M || 2 

E <*P (-3 2 r(d Xk h u - dy k h v ))\\ + Af 3 



1 - Us? 



M || 2 

M 1 - s v ) E CX P {-3 2 f(d x h u + d y K))\\ - Af ; 
\IU=i 

/II M || 2 

+ s„(s„ - s„) E ^p(-j^-d x h u )\\ - Af ; 
\llfc=i II 

/|| M || 2 

+(1 - s v )(s v - Su) E exp(i2-d Hfc ^) -AT 



— s^)(s t , — s u )- 



( ™ ( . ^d y h v \ M / . 27^ l. u \ « / 27r(d x h u + d y h v )\ 

exp ^ f ) }^ cxp ^ ) }^ cxp I * I 

27v(d x , h u +d y , h v ) 



M 




+ E 


exp 


fc = l 






M 


—M 


E 




k=l 




M 


-M 


E 




/ 



— ) E ox p (- 



) E ex P j 

k = l 



1 — UsNR 



s v (l-s u ) E exp (-j2f(dx fc /iu + dv fc 



Af- 



+Sv(s u - s v ) ( | E 0X P (-3^-d, Xk h u ) 
I fc— l 



M- 



+(l-s»)(s u ( | E exp(-j2fd I)fc h 1 ,) 

I fc — 1 



Af 



/ 



/ £ cxp (.-^^ff cxp fj^b^L) f cxp L. ^^ ) 

fc = l v ' k=l v 7 fc=l V / 

M / .2^d H h„\M / .2 w da. fc u \M { .2-*(d x h u +dy h v ) 
+ E exp ( —J f L cxp —J f I E exp J ^ S 

fc=i v ' k=i v y fc=i V 



V V 



-Af 
-Af 



^E exp ( - j h v ) 1 1 -« E^HHtM 



E exp (-jS 
fc=l 



. h u + d Wfc /i„)) +Af 3 



/ M 

1 + S U (1 - S U )U S NR I Af 2 — E CX P (-7^d xt; / ill ) 



1 + s„(l - s v )U s » 



E exp(-j^d„ 
Ifc=i 



(43) 



November 29, 2012 



DRAFT 



16 



and, of course, {G} Mt , = {G} vu . Consequently, the unconditional Weiss-Weinstein bound is 2 x 2 matrix given by: 
UWWB = HG _1 H T 

1 



[ 2 

I U V 



{G} 



(44) 



{G}uu{G}w - {G}( 

which has to be optimized over s u ,s v ,h u , and h v . Concerning the optimization over s u and s v , several other 
works in the literature have suggested to simply use s u = s v — 1/2. Most of the time, numerical simulations 
of this simplified bound compared with the bound obtained after optimization over s u and s v leads to the same 
results while their is no formal proof of this fact (see [9] page 41 footnote 17). Note that, thanks to the expressions 
obtained in the next Section concerning the linear array, we will be able to prove that s = 1/2 is a (maybe not 
unique) correct choice for any linear array. In the case of the planar array treated in this Section, we will only 
check this property by simulation. 

In the particular case where s u — s v — 1/2 one obtains the following simplified expressions 



2 1-^ -2(1-1^1) 1 



4 



M 2 



{G} uu — 



M 

E ex P {-3^d Xk h u ) 

k=l 



(l-%L 



USNR | 



M 



E exp(-j^d Xk h u ) 



2(1-%1)-2(1-|M(1 



USNR 

4 



k=l 



M 2 



-IT 



(45) 



{G} vv 



M 



E exp(-j^-d Vk h v 
fc=i 



(l - ¥ y 



Usmr M 2 _ 



M 

E e x P(-3T-d yk h v ) 

k=l 



-2T 



(46) 



and 



2 1 



M 2 — 



A I 



E exp {-j^{d Xk h u - d yk h v )) 



k=l 



-2 1 + M 2 



A I 



E ex P {-j^r{ d x k K + d Vk h v )) 
fe=i 



-T 



/ 



1 + M 2 



A I 



E ex P (-]T~d Xk h u ) 



k=l 



1 + M 2 



A I 



-T ' 



E exp(-j^-d yk h v ) 



k=l 



Again, the Weiss-Weinstein bound is obtained by using the above expressions in Eqn. d44l > and after an opti- 
mization over the test points. The optimization over the test points can be done over a search grid or by using the 
ambiguity diagram of the array in order to reduce significantly the computational cost (see lfl"8l . Il27l . [34|, |45|, 

ED). 
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2) Conditional observation model M.^: Under M2, let us set Csnr = -jr Y. 

" t= 

sions of the elements of matrix G are given by (see Appendix [F] for the proof): 

(l - ih £ 1 ) ex P ^s u (2s u - 1)Csnr [M - jr cos (^d Xk h u 

—d h 
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(48) 



(49) 



Af \ \ / / Af x 

exp ( 2s„(s M - l)CsAr fl ( A/ - J] cos {^fd Xk h u ) j j exp ( 2s^(s„ - 1)Csnr [M- cos {^f d Vk h v) 

(50) 



and {G} u „ = {G}„ u . Consequently, the conditional Weiss-Weinstein bound is 2 x 2 matrix given by using the 
above equations in Eqn. (l44l . As for the unconditional case, if we set s u = s v = 1/2, one obtains the following 
simplified expressions 
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2 (l - i ¥) - 2(1 - cxp ^-^p - 



{G} TO = 
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(53) 
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E cos (x^A) ~ E C0S (x d 2>A 

V V fc=l k=\ J J 

By using the above expressions in Eqn. (PHt and after an optimization over the test points, one obtains the 
Weiss-Weinstein bound. 



C. Source localization with a non-uniform linear array 

We now briefly consider the DOA estimation of a single narrow band source in the far area by using a non- 
uniform linear array antenna. Without loss of generality, let us assume that the linear array antenna lays on the 
Ox axis of the coordinate system (see Fig. [TJ, consequently, d Vi — 0, Vi. The sensor positions vector is denoted 
[d Xl . . . d XM ] . By letting 9 = sirup, where (p denotes the elevation angle of the source, the steering vector is then 
given by 



a(0) 



.2tt 

T 



dr,J 



.2vr , 
exp ( j—d XM t 



1 T 



(54) 



We assume that the parameter 6 follows a uniform distribution over [—1, 1]. As in Section HV-BI and since the 
parameter of interest is a scalar, matrix H of the test points becomes a scalar denoted hg. In the same way, 
there is only one element G [0, 1] which will be simply denoted s. The closed-form expressions given here are 
straightforwardly obtained from the aforementioned results on the planar array about the element denoted {G} UIt . 
We will continue to use the previously introduced notations Usnr — a -i (Ma\+a' } ) anc * ^snr = ^E I I s (Oil ■ 

1) Unconditional observation model M.%: The closed-form expression of the unconditional Weiss-Weinstein 



bound, denoted UWWB, is given by 
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A I 



M 

E ex P (-ij4 fc /ifi 
fe=i 

. \ -T 



-T 



E exp(-j^-d Xk hg 

k=l 



(55) 
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In order to find one optimal value of s that maximizes HG 1 H T , \fhg we have considered the derivative of 
HG _1 H T w.r.t. s. The calculation (not reported here) is straightforward and it is easy to see that aHG H 



ds 



0. Consequently, the Weiss-Weinstein bound has just to be optimized over hg and is simplified leading to 



(i-¥) ; 



USNR 

4 



Af 2 



UWWB = SUp- 
he 



M 

E ex P 1 

fc=i 



J A ax k 



he) 



-IT 



(56) 



2(1-%1) -2(1 



Usmr M 2 



M 

E ex P 1 

fc=i 



-j^rd Xk h g 



In the classical case of a uniform linear array (i.e., d Xk — d), this expression can be still simplified by noticing 

M 

that E exp (-j^d Xk h e ) =M exp (-j^h e ) . 

k=l 

2) Conditional observation model M.^'- The closed-form expression of the conditional Weiss-Weinstein bound 



CWWB is given by 



CWWB = 



h* ( 1 - M ) exp ( 4s(s - 1)C; 



M 



E cos (^-d Xk h e 

k=l 



I fee 



V 



exp 4s(2s - 1)C, 



M • 



M 

E 

fe=i 



E cos {^d Xk h e 



\ \ 



(57) 



exp 4(2s - l)(s - 1)C, 



M- E cos(^d ;Cfc /i e 



V 



-2(l-|ftfl|)exp 2s(s-l)C, 



M 



M 

E COS i 

k=l 



k=l 

-fd Xk he 



9HG H 



Again, it is easy to check that ^ 

HG _1 H T , Vhg is s = h. The Weiss-Weinstein bound is then simplified as follows 



0. Consequently, one optimal value of s that maximizes 



CWWB = sup 



\he\ 
2 



M 



exp ( -C SNR ( M - E cos (^-d Xk h g 



ho 2 (l-M) _ 2 (l-|^|)exp^-|C g ^ - cos (^-d Xk h 9 ) 



(58) 



In the classical case of a uniform linear array (i.e., d Xk — d), this expression can be still simplified by noticing 

M 

that E cos (^d Xk h g ) = M COS (^hg). 

k=l 



VI. Simulation results and analysis 

As an illustration of the previously derived results, we first consider the scenario proposed in [23 1 Fig. 5, i.e., 
the DOA estimation under the unconditional model using an uniform circular array consisting of M = 16 sensors 
with a half-wavelength inter-sensors spacing. The numbers of snapshots is T = 100. Since the array is symmetric, 
the performance estimation concerning the parameters u and v are the same, this is why only the performance with 
respect to the parameters u is given in Fig. [2] The Weiss-Weinstein bound is computed using Eqn. d45b . d46l > and 
(I47t . The Ziv-Zakai bound is computed using Eqn. (24) in |23l . The empirical global mean square error (MSE) 
of the maximum a posteriori (MAP) estimator is obtained over 2000 Monte Carlo trials. As in |23| Fig. (lb), one 
observes that both the Weiss-Weinstein bound and the Ziv-Zakai bound are tight w.r.t. the MSE of the MAP and 
capture the SNR threshold. Note that, in [23 1 Fig. (lb), the Weiss-Weinstein bound was computed numerically only. 
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Fig. 2. Ziv-Zakai bound, Weiss-Weinstein bound and empirical MSE of the MAP estimator: unconditional case. 



To the best of our knowledge, their are no closed-form expressions of the Ziv-Zakai bound for the conditional 
model available in the literature. In this case, we consider 3D source localization using a V-shaped array. Indeed, 
it has been shown that this kind of array is able to outperform other classical planar arrays, more particularly the 
uniform circular array ||47| . This array is made from two branches of uniform linear arrays with 6 sensors located 
on each branches and one sensor located at the origin. We denote A the angle between these two branches. The 
sensors are equally spaced with a half-wavelength. The number of snapshots is T — 20. Fig. [3] shows the behavior 
of the Weiss-Weinstein bound with respect to the opening angle A. One can observe that when A varies, the 
estimation performance concerning the estimation of parameter u varies slightly. On the contrary, the estimation 
performance concerning the estimation of parameter v is strongly dependent on A. When A increases from 0° to 
90°, the Weiss-Weinstein bound of v decreases, as well as the SNR threshold. Fig. [3] also shows that A = 90° is 
the optimal value, which is different with the optimal value A = 53.13° in |47| since the assumptions concerning 
the source signal are not the same. 

VII. Conclusion 

In this paper, the Weiss-Weinstein bound on the mean square error has been studied in the array processing 
context. In order to analyze the unconditional and conditional signal source models, the structure of the bound has 
been detailed for both Gaussian observation models with parameterized mean or parameterized covariance matrix. 

Appendix 

A. Closed-form expression of rj e (a, (3, u, v) under the Gaussian observation model with parameterized covariance 
Since y (i) ~ CM (0, R y (9)), one has, 

\K y (9)\ T ^ +f3 -^ 

)\ Ta \R v (6 + v)\ Tfi 

Q 



tt mt |R y (0 + u)| Tq |R y (0 



-^yp / ex p (-Ey H W rl y(')j rfY > ^ 
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Fig. 3. Weiss-Weinstein bounds of the V-shaped array w.r.t. the opening angle A. 



where r 1 = aK^iO + u) + f3R y l {8 + v) - (a + /3 - l)R y l {0). Then, since 



/ exp|-f:y^(*)r- 1 y(i)|^ = 7 r^|r' T 



one has 



Ve KA",v) 



\BL y (0)\ na+fl ~ 1) F\ T 



\R y (8)\ T(a+ ^ 



(60) 



\n y (e + u)\ la \R y (e + v)\ ip |R y (0 + u)r|R y (0 + v)nr-i| J 

B. Closed-form expression of rj g (a,/3, u,v) under the Gaussian observation model with parameterized mean 
Since y (t) ~ CAf (f t (0) , R y ) , one has 



(61) 



witlj] 



rj (a, (3, u, v) = — 1 j exp j (t) j dY, 



(62) 



£ (f) = a (y - f t (8 + u)f R y 1 (y - f t (8 + u)) + (3(y-f t (8 + v))" R y 1 (y - f t (8 + v)) 
+ (1 - a - p) (y - f t (8)f R y 1 (y - f t (8)) 
= y H R y 1 y+ai t H (8 + u) R^f t (8 + u) +pf* (8 + v) R y H t (6 + v) + (1 - a - 0) ff (8) R y % (8) 
-2 Re {y H R' 1 (af t {8 + u) + j3f t (8 + v) + (1 - a - 0) f t (0))} . (63) 



^ T5-1 



For simplicity, the dependance on t of f and y is not emphasized. 
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Let us set x = y - (af t (0 + u) + /3( t (0 + v) + (1 - a - ff) i t {6)). Consequently, 
x ff R y x x = y ff R y V - 2 Re {y H R y 1 (af ( (0 + u) + fSf t (0 + v) + (1 - a - f3) f t (0))} 
+ (aff (0 + u) + /3ff (0 + v) + (1 - a - 0) ff (0)) R y 1 (af t (0 + u) + /3f f (0 + v) + (1 - a - 0) f* (6) $64) 
And £ (i) can be rewritten as 

= x fl R y 1 x + (((), (65) 

where 

| (t) = a (1 - a) ff (0 + u) R y !f t (0 + u) +/3 (1 - /3) ff (0 + v) R y x f 4 (0 + v) 

+ (1 - a - /3) (a + fS) ff (0) R y x f t (0) - 2 Re {a/3ff (0 + u) R y J f t (0 + v) 

+ a (1 - a - 0) ff (0 + u) R y x f t (0) + (1 - a - 0) ff (0 + v) Rjff* (0)} . (66) 

Note that £ (t) is independent of x. By defining X = [x (1) , x (2) , . . . , x (T)], the function fjg (a, 0, u, v) 
becomes 

r) (a, (3, u, v) = ^ /t | r |T jf exp (~E X " R y ' x + * dX = ex P (*)J ' (67) 



since ^fer In exp ( - Ex H Ry : x ) dX =1. 



C. Closed-form expressions of |miR y 1 (0i) + TO2R y 1 (02)| end |TOiR y 1 (0i) + m2R y 1 (02) + J^Ry 1 (03)| 

Note that this calculation is actually an extension of the result obtained in |27| Appendix A in which mi = 
m 2 = \ and m 3 = 0, but follows the same method. The inverse of R y can be deduced from the Woodbury formula 

1 / ^a(0)a"(0) \ 



a2||a(0)|| 2 +< 
Then, 



Since the rank of &(9k)& H (0fe) is equal to 1 and since 0i 7^ 02 7^ 03 (except for = h; = 0), the above 

3 

matrix has M — 3 eigenvalues equal to \ ^2 m k and 3 eigenvalues corresponding to the eigenvectors made from 

CT "fe=i 

the linear combination of a(0i), a(02), and a(03): a(0i) +pa(02) + qa(03). The determinant will then be the 



product of these M eigenvalues^. Let us set 



^||a(0 fc )|r + <r 

Then, the three aforementioned eigenvalues denoted A must satisfy: 



Vk = „ ,„ :„ 2 , _ 2 . k = 1,2,3. (69) 



J]m fe Ry 1 (0 fc ) (a(0i) + pa(0 2 ) + ?a(0 3 )) - A (a(fli) + pa(0 2 ) + 9 a(0 3 )) • (70) 



\k=l 



4 Note that we are only interested by the eigenvalues. Consequently, the lineal" combination of of a(0i), a(<?2), an d a (0:i) can De written 
a(0i) +pa{0 2 ) + <?a(0 3 ) instead of ra(fli) + pa(0 2 ) + ?a(0 3 ) 
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By using Eqn. d68l > in the above equation and after a factorization with respect to a(<?i), 9.(82), an d a (f? 3 ) one 
obtains 

(x - m lVl \\a(9 1 )\\ 2 -pm lVl a H {9 1 )a{9 2 ) - qm Wl a H (0i)a(0 3 )) a(9 1 ) 

-m 2V2 a H (e 2 )a{e 1 ) +p(x- m 2 ^ 2 ||a(0 2 )|| 2 ) - qm 2 ^ 2 a H (0 2 )a(0 3 )) a(0 2 ) (71) 
^ -m 3 ip 3 a H (9 3 )a(9 1 ) - m 3 <p 3 pa H (8 3 )a(0 2 ) + q(x- m 3 ^ 3 ||a(0 3 )|| 2 )) a(0 3 ) = 0, 

wherjf] 

x = 1 - <T 2 n \. (72) 

Consequently, the coefficients of a(9i), a(9 2 ), and a(0 3 ) are equals to zero leading to a system of three equations 
with two unknown (p and q). Solving the two first equations to fincQ p and q, and applying the solution into the 

last equation, one obtains the following polynomial equation of x 

3 33 
x 3 _ JYjnWk \H^)f " f EE m *VW^fc< (||a H (0 fe )a(0 fe ,)|| 2 - ||a(0 fc )|| 2 ||a(0 fc ,)f) (75) 

k=l k=lk'=l 

k'^k 

-mimamaMtipa (||a(0x)|| 2 ||a(0 2 )|| 2 ||a(0 3 )|| 2 - ||a" (0 2 )a(0 3 )|| 2 \\a(0x)f (76) 
- \\a H (8 1 )a(e 2 )\\ 2 ||a(0 3 )|| 2 - \\a H (9 3 )a(9 1 )\\ 2 \\a H (9 2 )\\ 2 + a H (9 3 )a(9 2 )a H (0 1 )a(0 3 )a H (0 2 )a(0 1 )(77) 
+a if (0 3 )a(0 1 )a" (0^(0^ (0 2 )a(0 3 )) = 

Since we are only interested by the product of the three eigenvalues, we do not have to solve this polynomial 

3 

in A and only the opposite of the last term is required. This leads to Eqn. OTJ with m k = 1. Of course, the 

fe=l 

2 

closed-form expression of |miRy 1 (0i) +m 2 R~ 1 (0 2 )| is obtained by letting m 3 — and = 1 in Eqn. 

fc=i 



D. Closed-form expressions of £ e (/Li, p) 

Remind that the function C, e (fi, p) is defined by Eqn. ( TT8T >. Let us define p as the number of parameters per 



sources (assumed to be constant for each sources). Then, without loss of generality, the full parameter vector 9 

r 1 T 

can be decomposed as 9 — 9\ . . . 9 T N where 9i = [O^i . . . 6i :P ] , i — 1, . . . , N with q = Np. Remind that 

T r 1 T 

p = [0 ... /ij ... 0] and p = [0 . . . Pj . . . Oj .It exists two distinct cases to study: when both index i and j are 



3 

5 Note that, from Eqn. I16K mj, = 1. 

k-l 

6 p and g are given by 



m 2 y? 2 a w (0 2 ) 


(m m a(0 1 )a H (0 1 ) + 


— mnp 1 ||a(fli) 


| 2 )l)a(0 3 ) 




(m 2¥ > 2 a(0 2 )a"(6> 2 )4 


(x — m 2 (/3 2 |a(0 2 ) 


| 2 )l)a(0 3 ) 



and 

(x - m^j |a(0i)|| 2 ) (x - m 2 ip 2 ||a(0 2 )|| 2 ) - m 1 tp 1 m2<p 2 a - H ( i) a ( s 2)a H (0 2 )a(0i) 



mi Vl a»(0i) (m 2V2 a(0 2 )a»(0 2 ) + (as - m 2 <^ 2 ||a(0 2 )|| 2 ) i) a(0 3 ) 



(73) 



(74) 
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such that (m — l)p + 1 < i < mp, m = 1, . . . , N and (m — l)p + 1 < j < mp or when (m — l)p + 1 < i < mp, 
rn = 1, . . . , N and (n — l)p + 1 < j < np, n = 1, . . . , N with m ^ n. Therefore let us denote: 



p m = [o---o hi o---ofer 
Pm = [o---o hi o-ofer 



if (m — l)p + 1 < i,j < mp 



and 



p m = [0---0 /ii 0---0] T eW, 
p„ = [0---0 h 3 0---0] T eK p , 



if 



(m — l)p + 1 < i < mp 
(n — 1) p + 1 < j < 



, with m 7^ n. 



7) 77;e case where (m — l)p + 1 < i, j < mp: In this case, one has: 

A(0 + p)-A(0 + p) = [O---O a(0 ro +pj-a(0 ro +p m ) • • • 0] e C pxA \ 
and consequently, 

T 

C e (A*, P) = R» 1/2 (a (0 m +M m ) - a (0 m +pJ)| £ ||{s (i)} J| 2 . 



Due to Eqn. d28l ). one has 

R n 1/2 (a(0 m +pj - a(0 m +p m )) 



M M 



EE { R n 1 }„ ex p ( j't ^ t - r H 0m ) ( exp ( -^T r ^ 

i=lj=l 



27T 



2tt 



2tt 



x ( exp ( jyr^ p m j - exp I j— Tj p m 



In particular, in the case where R n = cr^I one obtains 

2 1 M 

R n 1/2 (a(0 m +/x m ) -a(0 m +pj) = 



2tt 



exp 1 -J— r i Mm I ~ exp I -j— v % p m 



2tt 



(78) 



(79) 



(80) 



(81) 



(82) 



(83) 



2) The case where (m ~ l)p + 1 < i < mp and where (n — l)p + 1 < j < np: Without loss generality, we 
assume that n > m. Then, 

A(0 + p) - A(0 + p) = [a(0i)-a(0i)---a(0 rri + /x m )-a(0 m )---a(0„)-a(0 n +p„)---a(0 2V )-a(0 iV )] 
= [0 a(0 m +pj-a(0 m ) a (0 ro ) - a (0„+p n ) 0-0], (84) 

and consequently, 



C e (p, p) = E || R n V2 ( a ( ™+P m ) - a (0 ro )) {s (i)} m + (a (0 n ) - a (0„+pJ) {s (*)}, 



(85) 



November 29, 2012 



DRAFT 



25 



Let us set x. = R n 1/2 (a (0 m +/Lt m ) - a(0 m ))and £ = R„ 1,z (a(0 n ) - a(6 n +p n )) . Then, Co (M>p) can be 
rewritten 

T 

C e (/,,p) = ^||^{s(t)} m + e {s(i)}J| 2 
t=l 

= £ (>* ||{s (i)}J| 2 + x h q{s {t)}* m {s (t)} n + e H x{ s (t)} m {s (t)>; + ||{ s (*)}„ 



_ u-i/a. 



= ||{s W} m f + Hi 8 Win I! 2 + 2 Re >c H e J2 (*)>; {s (t)}„ • (86) 

By using the structure of the steering matrix A, it leads to 

M M 

* H k=EE {V},, exp (if (rj - rf ) m ) exp (-if rf Mm ) exp (if rj/x m ) , 
^=EE {Rn 1 }, , exp (if (rj - if) n ) exp (-if rfp n ) exp (if rjp n ) , (87) 

i=lj=l 
M M 

A^EE {Rn 1 }, , exp (if (rje n - vj6 m )) exp (-if rf/xj exp (if rjp n ) . 

i=ii=i 



£. Proof of Eqn. (HD, and dH 

In fact, one only has to prove Eqn. ( l43l l since Eqn. fiTt and ( |42T > can be obtained by letting = and s u = s„ 
in Eqn. d43l l and by using s u ) for Eqn. ( HTb and (h v , s v ) for Eqn. (1421 1. By plugging Eqn. ( f30b and d32l i into 
Eqn. ( fTSI ), and by considering the following expressions 



M 



a H (8 + h u )a(8 + h„) = J2 exp (if (<2 y ,/i„ - d Xz h u )) = {a H (8 + h v )a(8 + h u )) 



M 



M 



a H (8 ± K)a(8) = £ exp (T3%d Xi h u ), and a ff (9 + h u )a(9 - h u ) = £ exp (-jfd Xi h u ) 

i=l i=l 

one obtains the closed-form expressions for the set of functions i) g (a, f3, u, v) 



Ve{ s uj s »j h u , h u ) — 



1 — USNR 



S u Si- 



A I 



E ex P {-3^{d Xk h u - d Vk h v )) 



k=l 



_ t - 5 1i (l S u SyJ 



+s v (l - s u - s v ) 



T 2 2 



M 

E exp (-if d Xfc /i„ 
fe=i 

A/ 

E exp (-jf 

fc=i 



- M 2 

2 ' 



M 2 



M 2 



U a' 
M 



E( . 2nd y h v \ M / .2nd m h u \ M ( . 

exp I j — f — E ex P ( -J — a — E exp I J 

k=l v 7 fc=l v 7 fc=l v 

E' / ■ 2ixd y h v \ M ( .2ird w .h u \ M 

exp ( -j — f— I 2^ exp ( j — f — I L exp 

fe=l v 7 fe=l v 7 k=l 



k= 

-M 
-M 



27r(d Xk h u —dy k hv) 



. 2ir(dx k h u —dy k h v ) 
'3 A 



M 



E exp (-if dy k h v ) 



k=l 
M 



M 



E exp (-if (d Xk h u - d Vk h v )) 

k=l 



E exp (-if d Xk h u ) 
M 3 



fc=i 

2 



(88) 
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/ / 



1 - UsNR 



(1 - s u )(l - s„) 



M 



E ex P 0^r(dx k h u - d Vk h v )) 



k=i 



+(1 - s u )(s u + s v - 1) i E exp (j^d Xfe /i„) 



+{1 - s v )(s u + s v - 1) I J2 ex P(j 2 T d yk h v ) 




\ \ 

— M 2 
M 2 ) 



-(1 - «„)(! - s v )(s u + s v - l)^%^x 



V 



2lT(da: k h u -dy k h tl ) 

A 



E exp ( /^M £ exp (-j 27r ^ fcfe " ) E ex P (j 

k=l v y k=l v ' fe=l v 

^i, / .27rd wt h«\ ^ / . 2nd x . h u \ / . 2n(d Xl h u — d y , h v ) \ 

+ E cx p (-j — 5 s — J E ex p(j — ) E cx p(-j A 
fe=i v ' fe=i v y fe=i v y 



-M 
— M 



M 



E exp(-j^d !/)b /i„) 



fe=i 

M 



M 



M 



E ex P {-i^{d Xk h u - d yk K)) 
k=l 



E ex P \-j-^d Xk K 

+ M 3 



fc=l 

2 



/ / 



1 - UsNR 



+ (1 - s„)(s„ - s„) 



M 

E ex P {-j^f{dx k K + d yk h v )) 
fe=l 

2 

E exp(-jf d^fcj -M 2 
fe=i 

M 

2tt 



— M 2 



E expOx^M 



fc=i 



-M 2 



E" / .2wd y ,h v \ ^ ( . 2nd Xk h u \ J^, ( 
exp I J cxp M j*— E ex P 

fe=l v 7 fe=l v y fe=l v 

/ . 2nd y ,h v \ ^ / . 2nd Xk h u \ j^, ( . 

+ E cx p ( — f— E cx p ( — t — ) E ex P [J 

k=l v 7 fe=l v ' k=\ v 



. 2ir(d Xk h u -\-d t 



27r(d !Cfe h u +dy k h v 
A 



-M 
-M 



M 

E e x V>(jird yk h v 
fe=i 

M 

2_7r 
A 



M 



M 



E exp (-j^<i Xfe /i, u ) 



fc=i 



E ex P {-jjr(d Xk h u + d yk h v )) 
fe=i 



+ M 3 
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%(1 s u ,s v , h u ,ht,) 



and 



1 — UsNR 



V 



~\~S V (S U Sy J 

+(1 - s u )(s u - s„) 



M 

E ex P (-jjc( d x k K + d Vk h v )) 
fe=l 

2 

fe=i 



— M 2 I * 



/ M / , \ M 



M 

E ex P (-j^d y)t /i„) 

k=l 



E( .2irdy.hv \ £\ ( ■ l-itdx Ji„ \ My ( . 2n(d as h u +d y h v 
exp ij — f— }J exp I j — ^ £ exp I -j ^ 

k=l v y fe=l v ' k=l v 



M 2 



exp ( — j — ) exp ( — j 



M 

+ E 

/c=i 

— M 



E 

k=l 



.2-Kd x ,h u \ / .27r(d x , h u ~\-d v ,h v ) 

J2 exp I j ^ s — 

k=l v 



-M 



M 



M 

E exp (-j 221 

fe=i 

E ex P {-i^rid^K + d yk h v )) 
fe=i 



a/ 



E exp(-j 2 ^d a;fc / l „ 
M 3 



fe=i 

2 



%(s«, 0, h u , 0) = 1 + s„(l - s u )i7 ; 



fc=i 



2tt 



?7 e (0, s„, 0, h„) = 1 + s v (l - s v )Usnr M 2 



— T 



-2' 



(91) 



(92) 



(93) 



One notices that the set of functions fj g (a, (3, u, v) does not depend on 6. Consequently, it is also easy to 
obtain the Weiss-Weinstein bound (throughout the set of functions r\ (a, (3, u, v)) by using the results of Section 
IIV-BI whatever the considered prior on 9 (only the integral f & ^+J^rj^jd8 has to be calculated or computed 
numerically). In our case of a uniform prior, the results are straightforward and leads to Eqn. ( |4TI ). ( l42l and ( l43l ). 



F. Proof of Eqn. gU), g9|) and (O 

The set of functions -q e (a, f3, u, v) is given by Eqn. (I17t . So, it only remains the calculation of functions £ e (/it, p) 
from Eqn. (1181 . Since R n = er^I, one obtains 
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C e (h n ,0) = Ce(-h»,0) 
C e (h„,0) = Co(-h»,0) 
C e (hu)-h«) = Ce(-h«,h u ) 
Ce(ht>i-h«) = Ce(-h K ,h„) 

< 

CeChu,!^) = Ce(ht,,h u ) 

= Ce(-h„, -h„) 
C e (-h„,h„) = Ce(h„,-h.„) 
= (e(- h v,K) 
( e (h u ,h u ) = Ce(h„,h„) 
Again, since the set of functions C, g (fi, p) does not depend on 6, the set of functions i] g (a, (3, u, v) is given by 
plugging the above equations into Eqn. (fTTT i and does not depend on 9. Consequently, as in unconditional case, the 
set of functions rj (a, f3, u, v) is obtained by using the results of Section HV-BI whatever the considered prior on 9. 
In our case of a uniform prior, the results are straightforward and leads to Eqn. (1481 . d49b and (1501 , 
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